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We describe a class of exactly soluble models for gravitational collapse in spherical symmetry 
obtained by patching dynamical spherically symmetric exterior spacetimes with cosmological interior 
spacetimes. These are generalizations of the Oppenheimer- Snyder type models to include classical 
and quantum scalar fields as sources for the interior metric, and null fluids with pressure as sources 
for the exterior metric. In addition to dynamical exteriors, the models exhibit other novel features 
such as evaporating horizons and singularity avoidance without quantum gravity. 



INTRODUCTION 

Gravitational collapse of matter configurations has 
been an important area of investigation for many years 
due mainly to the interest in understanding the dynam- 
ics of black hole formation. There are only a handful 
of analytical models in general relativity beginning with 
the very first by Oppenheimer and Schneider (OS) [ij. 
In their model stellar matter is assumed to be homo- 
geneous and isotropic and modelled by a closed FRW 
solution with pressureless dust. This was followed some 
years later by the Vaidya solution [2| which describes the 
collapse of pressureless null dust to a Schwarzschild black 
hole. Models of this type have been much studied in var- 
ious contexts [3Hl] 

There is a generalization of the Vaidya solution to a 
null dust with pressure, with equation of state P = kp 
[8|, which gives the Reissner-Nordstrum charged black 
hole for /c = 1, and a more general class of hairy black 
holes for k > 1 as the end point of collapse. Although 
these solutions are analytic, they describe purely ingoing 
(or outgoing) matter in spherical symmetry. 

The simplest asymptotically flat solutions with both 
inflow and outflow in spherical symmetry were found nu- 
merically, with a minimally coupled scalar field as the 
matter source [9 . These solutions exhibit a rich struc- 
ture, including a discrete self-similarity and a mass scal- 
ing law for black holes. This work has led to much ad- 
ditional research on gravitational collapse [10 , including 
quantum gravity inspired effects at the onset of black 
hole formation llTHMl . The latter are expected to play 
a fundamental role in the late stages of collapse and it 
is not unreasonable to expect that quantum gravity will 
drastically affect the strong field regime. Indeed these 
works show that black holes form with a mass gap, a re- 
sult that was also suggested in an OS type model with 
a quantized scale factor in the interior region [15]. The 
quantum gravity corrections used in the numerical sim- 
ulations of scalar field collapse were inspired by the so- 
called polymer quantization procedure [16^ which grew 
out of the loop quantum gravity. 

An interesting feature of polymer quantization is that 
it introduces a length scale in addition to h into the quan- 



tum theory. This comes about because the choice of 
Hilbert space used to represent operators is such that 
conventional momentum operators (ie. generators of 
translations) do not exist, but are defined indirectly as 
certain functions of translation operators. This is be- 
cause translation operators are not continuous in the 
translation parameter, unlike in Schrodinger/Heisenberg 
quantization. A direct consequence is that kinetic en- 
ergy operators are bounded above. This may be viewed 
as an ultraviolet cutoff that is built in due to the choice of 
Hilbert space. It is this fact that has the potential to lead 
to interesting new physics at short distances while recov- 
ering standard quantum mechanics at large distances. 
This quantization prescription has been applied to the 
scalar field [17-20 and to semiclassical gravity in the 
case of homogeneous cosmology [21 . 

Motivated by these works we describe an analytical 
model of gravitational collapse that combines three fea- 
tures: (i) the basic OS idea, (ii) the generalization of the 
Vaidya metric [8 and (iii) polymer quantization of the 
matter sector. The last of these has been investigated 
in detail recently in the quantization of the scalar field 
on a Friedmann- Robertson- Walker (FRW) background, 
leading among other things, to the remarkable prediction 
of an extended inflationary period in the early universe 
without a mass term or other scalar potential [21 . Our 
motivation for using this quantization is partly motivated 
by this result. 

The model we explore is a classical FRW interior met- 
ric sourced with a polymer quantized scalar field that is 
patched using junction conditions to the generalization of 
the Vaidya metric given in [8 . The setting is thus that 
of a quantum field on a fixed classical background, and 
so does not incorporate any quantum gravity corrections 
(unlike Ref. [TF which has some of the same ingredients). 
Our main result is that polymer quantization of matter is 
sufficient to avoid a curvature singularity - surprisingly 
without recourse to quantum gravity. 

In the next section we describe the classical model. 
This is followed in section III by a discussion of the poly- 
mer quantization of matter and its effect on the dynam- 
ics of the FRW scale factor. In Section IV we apply this 
dynamics to study collapse scenarios using the junction 
conditions at the interface of the interior and exterior 
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metrics. The concluding section contains a discussion of 
the main results and its relation to some other works. 



THE MODEL 

The collapse models we consider are all constructed by 
patching together an interior FRW spacetime with a gen- 
eralized Vaidya-like exterior. This differs from the usual 
OS model in that the interior solution is to have quan- 
tized scalar field matter, and hence non-zero pressure. 
Therefore matching to a pressureless null dust is not pos- 
sible without surface stresses; this is why we must use a 
more general exterior solution. 

The parametrization we use for the interior metric 



-de 



aHt) 



{dr^ + r^d0^ + r^sm^ed<p'^), 

(1) 



{l + r2/4)2 

and that for exterior metric (/^^ is 

dsl = -f{v, R) dv^ + 2 dvdR + {d0^ + sin^ Odcp^) , 

(2) 

which is written in advanced Eddington-Finkelstein coor- 
dinates (v^r). These coordinates are convenient because 
the trapping horizon, which is one of the objects of inter- 
est here, is given simply by f{v,R) = 0. We have taken 
the coordinates {0, (j)) on the 2— spheres to be the same 
in the two metrics. The coordinates (t, r) of the interior 
metric and (v, R) of the exterior metric are of course dif- 
ferent and so the matching of these metrics along a com- 
mon timelike 3— surface X) must take this into account. 
For this we need the induced metrics h^^ and extrinsic 
curvatures K^^ of E from both sides to carry out the 
standard matching analysis. 



Matching surface: interior view 

Let the metric on E from the FRW side be given by 
setting r = ro, a constant. This is the natural choice that 
describes a 2— sphere evolving along a timelike trajectory. 
The metric is then 



(l + roV4)2 



((i6>^ -h sin^ 6>#^) . (3) 



The unit timelike tangent and unit spacelike normal of 
the surface H are respectively 



dt 



a \ A I \or 



(4) 



Using these, the non-zero components of the extrinsic 
curvature 



(5) 



are 



K2 



l-r|/4^ 
a{t) ro 



(6) 



Matching surface: exterior view 

Let the timelike surface from the exterior side be given 
by i? = R{t) and v = v(t), where t is the interior's proper 
time coordinate. The induced metric /i^;, on S is then 

dsl = h+^dx'^dx'' 

= - {f{v{t), R{t))v^ - 2i?v) df 
+R{tf {d9^ + sin^ edcj)'^) , 



(7) 



where the dot denotes d/dt. 

The unit timelike tangent and unit spacelike normal of 
H are respectively 



.a 



1 



/i;2 - 2Rv L 



R 



d_ 

M 



(8) 



^fv'^ -2Rv 
The extrinsic curvature components 



R{dv)a^v{dR)a . (9) 



(10) 



are therefore 



^tt — ^ah^+^+ 



(flRV - Ivv - ^IrR) + ^Rv - Rv) 



2 (^fv'^ - 2Ri?j 
where f^R = df/dR etc. and 



. \ 3/2 



fv-R 



(11) 



(12) 



R\ fv^ - 2Rv 



Junction conditions 

The Israel junction conditions require continuity of the 
metric, and continuity of the extrinsic curvature if there 
is to be no surface stress-energy. Taken together these 
equations describe the dynamics of the boundary surface. 
Matching the induced metric components gives 



a{t)b = R(t), 



2Rv = 1, 



(13) 
(14) 
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where we have set b = ro/(l + ro/4), which is a constant. 
Matching the extrinsic curvature components, and using 
the previous equation give 

fv-R = c, (15) 

{{ff,R - f,v)v - 31rR) 

+ 2{Rv -Rv) = (16) 

where c = (1 - rg/4)/(l + r^/4); note that 6^ + = 1. 
The last condition may be rewritten in a simpler form by 
noting that the derivatives of and (jTsl) give 



2Rv = 
which leads to 



2vR = fv{2R- fv), (17) 



2R 



0. 



(18) 



The four equations ( T3p5 ) and (18) fully determine the 
dynamics of the boundary and the exterior metric func- 
tion /(i?, v): a matter source for the interior FRW deter- 
mines a(t), and hence R{t) through eqn. (13), the next 



two determine f{R^c) and v{R^c)^ and the last equation 
requires / = f{R) on the boundary. The second and 
third equations also give the four- velocity of the bound- 
ary as seen from the exterior, ie. 

:= (i;,i^,0,0) 

= 0(c±y^^),±y^^,O,Oy (19) 

The ± solutions correspond to collapsing and expanding 
solutions. We note also that = implies that are 
tangent to radial geodesies in the exterior spacetime. 

We are interested in the function R{v) which gives the 
trajectory of the FRW boundary. From (19) this is 



dR 

dv 



7^ = ±V^f (c T V^) (20) 



This is equation we study for various cases of interest 
for the classical and quantum scalar field. Its use re- 
quires only the function f{R) obtained from the junction 
condition, which in turn depends on the interior matter 
Hamiltonian. If the interior metric is taken to be flat 
FRW the matching conditions above remain valid with 
the changes b = tq and c = 1. We restrict attention to 
this latter case for explicit calculations and comment on 
the other values c < 1 in the discussion section. 

We will see that for scalar field matter in the interior, 
the exterior metric function determined on the boundary, 
/(i?), has a dynamical extension which can be analyti- 
cally determined and interpreted as a null fluid with pres- 
sure. Furthermore, if the interior matter field is quan- 
tized using polymer quantization, which is a semiclassical 
approximation, the surface trajectory exhibits qualita- 
tively new features. Both classical and quantum matter 
cases exhibit interesting dynamical horizon behaviour. 



CLASSICAL MATTER SOLUTIONS 

Our goal is to construct models of gravitational col- 
lapse that go beyond the OS solution by using quantum 
matter in the interior. One way to do this is to begin 
with the Hamiltonian formulation of the gravity-matter 
dynamics in the interior and then quantize the matter 
canonical variables. This is a canonical approach to the 
semiclassical approximation. 

We begin with the Arnowitt-Deser-Misner (ADM) 
canonical form of the 3+1 action for Einstein gravity 
minimally coupled to a massless scalar field 



j dt j d^x [it' 



Qab ■ 



Nn-N^Ca), (21) 



where (tt^^, Qab) are the ADM canonically conjugate vari- 
ables and H and Ca are the Hamiltonian and diffeomor- 
phism constraints, and (^,^0) are the scalar field canon- 
ical variables. Reduction to homogeneous and isotropic 
cosmology is attained by the parametrization 



Qab = a^(t)eab, 



~ab ^ P^^ab^^ 

2a 



(22) 



where Cab is the flat spatial metric in Eqn. ([T]). This 
gives the reduced action 



S = Vn dt 



P 
"24a 



Stt 



(23) 



for flat FRW in the interior (which corresponds to c = 1 
and 6 = ro in the junction conditions), and Vq = J d^x 
comes from the reduction to homogeneity. 



Dust: Oppenheimer- Snyder solution 

From a Hamiltonian perspective this model may be 
viewed as arising from the dynamics of the canonically 
conjugate pair (a,p). Evolution is given by the Hamil- 
tonian constraint above but with the scalar field energy 
density term replaced by a constant: 



n = 



P 

24a 



167rpo = 0. 



(24) 



where po = constant. The canonical equations of motion 
for lapse = 1 are 



P 



{a,n} 

{p.n} 



p 

'l2a 
24^* 



(25) 
(26) 



Rewriting the constraint equation \2A\ using the d evolu- 
tion equation, the junction condition R{t) = roa{t) and 
the R equation (Il9| (with c = 1) gives 

W^po = 6ad' = ^(1 - f{R)), (27) 
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which gives the known result 



47r 



(28) 



which constrains the exterior metric to be Schwarzschild 
and provides the interpretation that the interior is dust 
bah of radius tq and density pQ. 



Scalar field 

As a warm up to the quantum matter problem, it is 
useful to consider the interior FRW sourced with a mass- 
less minimally coupled scalar field. Although this is a 
natural generalization, it doesn't appear to exist in the 
literature. As we will see the exterior metric takes an 
unusual form, but is one of a class of known solutions. 



Hamilton's equations obtained from the action (21) 
with N = 1 are 



a 

P 

P0 



P 

12a 

p2 24^FoP^ 



24a2 



WttVoPc/) 



together with the Hamiltonian constraint 



n 



P 
"24a 



0. 



(29) 

(30) 
(31) 
(32) 



(33) 



The procedure for finding the exterior metric function 
on the boundary uses the Hamiltonian constraint, the 
equation of motion for the scale factor, the junction con- 



dition R{t) = ba{t)^ and the trajectory equation (19). We 
have 



= n{p^,RJ{R))=0. 



(34) 



For the Hamiltonian constraint (33) this gives 

4 



SttVopI 



6 



R' 



'0 



6|5 



(35) 



This determines the exterior metric function on the 
boundary to be 



f(v,R)\, 



1 



3 i?4 



1 



a 



(36) 



where a > is a constant since is a constant of the 
motion. This suggests unusual exterior spacetimes as 
extensions of this boundary function. 



Exterior solution 

It is interesting that there is a large class of exact solu- 
tions with null fiuid stress energy tensor [8 that provide 
possible extensions to the exterior of the boundary func- 



tion (36). These are given by the metric 



ds' = 



2M{R,v) 
R 



dv^ + 2dvdR + R^dQ^ 



which arises from the source 
1 dM 



VaVb 



27ri?2 dv 

-\-p{v, R)W(^aVb) + P{V^ R) {gab + ^(a^6)) • 



(37) 

The pressure P{v,R) and energy density p{v,R) are given 

by 



4^^2/c+2 



where 



and 



M{R,v) m{v) 



2{2k - l)i?2/c- 



Va = (1,0,0,0), 
Wa = (F/2, -1,0,0), 



(38) 



(39) 



(40) 



are the future pointing null vectors; F is the coefficient 
of dv"^ in the metric, and /c is a real parameter. 



From the form of the metric (37) it is apparent that 



there is a unique value, /c = 2, that extends the boundary 
function f{R) = 1 — a/R^ derived in eqn. (36) to a static 
exterior. For this case we must have m(v) = and g{v) = 
constant. 

There are however interesting dynamical possibilities 
for the exterior for other values of k. These are obtained 
by finding the functions m{v) and g{v) such that on the 
interface R{v) given by ([2Q|) we have 



1 



2m{v) 

R(v) '^J2k- l)R{v) 



2k 



1 



R{vY 



and 



l{v) 



2m{v) 

R{v) ^ J2k- l)R{v) 



2k ' 



(41) 



(42) 



The first of these is the requirement that the exterior 
metric match the boundary function, and the second is 



the junction condition (18). Taking the derivative of the 



first and using the second gives the unique one parameter 
{k) family of solutions 



m{v) 



a{2k-A) 
{2k-l)R{vf'' 
-3aR{vf^-\ 



(43) 
(44) 



FIG. 1: A: = 1 exterior with classical scalar field interior: the exterior spacetime is dynamical and has null singularity at a fixed 
value of advanced time (here vq ~ 0.67). 




FIG. 2: k — 2 exterior with classical scalar field interior: the only solution is a static exterior (constant m{v) and g{v)). This 
is qualitatively similar to the Oppenheimer-Snyder solution. 



The case k = 2 gives m{v) = and g{v) = —3a, which is 
the static solution already noted above. All other values 
of k provide a dynamical exterior. 

Typical classical matter solutions are exhibited in the 
Figures 1-3. Each shows the junction trajectory, the met- 
ric functions m{v) and g{y) and the corresponding Pen- 
rose diagram with a (dynamical) horizon trajectory. All 
calculations are with the initial condition R(v = 0) = 2 
and the parameter values = 1 and ro = 1. 

k = 1 : (Fig. 1) A dynamical horizon forms and evolves 
on a time like trajectory and then becomes null and tan- 
gent to the surface trajectory; there is a curvature singu- 
larity at the corresponding value of v. The horizon starts 
to shrink because the metric functions are such that an 
ingoing flux of positive energy evolves to an ingoing flux 
of negative energy. This is apparent from examining the 



stress-energy tensor (37). The null singularity is naked 



because timelike observers in the exterior region intersect 
the V = Vq line in finite proper time. 

k = 2 : (Fig. 2) This is the static case and is qualitatively 
similar to the Oppenheimer-Snyder solution: m{v) — 
and g{y) = constant < which gives a single horizon. 
The Penrose diagram is therefore the standard one rep- 
resenting gravitational collapse to a black hole. 

k = 3 : (Fig. 3) This is another case with dynamical ex- 
terior and has the surprising feature that the dynamical 
horizon becomes null at finite v and extends to spatial 
infinity. The change from k = 1 to k = 3 dramatically 
changes the horizon dynamics, as may be expected since 
the order of the equation F{r^v) = is quite different. 
The metric functions are such that the null energy con- 
dition holds. 
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QUANTUM MATTER SOLUTIONS 

In this section we consider the same model as in the 
last section, but with a quantum scalar field. The re- 
sulting model is in the context of the usual semiclassical 
approximation where the expectation value of quantized 
matter is used as a source for the classical gravitational 
field. 

In this approach the interior dynamics is derived from 
the semiclassical constraint 



24a 







(45) 



where is suitably chosen state could depend on the 
scale factor a{t) and other parameters. A guiding princi- 
ple for selecting the state is that the resulting spacetime 
give the expected classical result for a large universe, but 
with increasing quantum affects as it shrinks in size and 
approaches the would be classical singularity. This is the 
type of state we select, and it does depend on the scale 
factor. 

The derivation of the boundary function is identical to 
that for the classical case since the junction conditions 
do not depend on the nature of the interior dynamics. 
The semiclassical constraint gives the formula 



f{v,R)\j: = l 



{^\pimR), (46) 



where the scale factor (and hence R) dependence of the 
expectation value is made explicit. As we will see, it is 
this dependence which can drastically modify the exterior 
solution. 

The quantization procedure we follow is that described 
in [21 , where a particular type of polymer quantization 
is applied to the scalar field. The variables used in this 



approach are motivated by but different from the ones 
that arise from loop quantum gravity [16 . The essential 
difference is the manner in which the field translation 
operator is defined. 

In the following we review this quantization procedure 
and summarize how it leads to a modified interior so- 
lution. The main feature of interest is that the interior 
scalar field energy density turns out to be bounded above; 
matching this new interior solution to the exterior solu- 



tion (37) gives an evaporating horizon and a remnant as 



the end point of gravitational collapse. This is one of the 
main results of this paper. 



Polymer quantization 

This quantization method does not use the standard 
canonical variables ((/), P^) as the starting point of quan- 
tization, but rather the new variables 



(t)f = I c^^xv^/(x)(/)(x), Ux exp 



J 



, (47) 



where the smearing function /(x) is a scalar [21]. The 
parameter A is a spacetime constant with dimensions of 
(mass)~^, and the factor in the exponent is required 
to balance the density weight of p^. These variables sat- 
isfy the Poisson algebra 



{4>f,Ux} = ifm. 



(48) 



Specializing to a scalar field on an FRW background, 
these variables become 



(/)/ = Voa^cj), Ux = exp (iXp^/a^) , 



(49) 



where Vb = J d^x is a fiducial comoving volume and / is 
set to unity since this is a reduction of (/)/ to the spatially 
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FIG. 4: /c = 2 exterior with quantum scalar field interior: the shell trajectory exponentially approaches r = and there is no 
curvature singularity. The apparent horizon forms and evaporates in finite time. 



homogeneous case. The Poisson bracket of the reduced 



variables is the same as that of the unreduced ones (48). 



Quantization proceeds by realizing the Poisson algebra 



(48) as a commutator algebra on a suitable Hilbert space; 
the choice for polymer quantization has the basis |/i G 
R} with inner product 



(50) 



where 5 is the generalization of the Kronecker delta to 
the real numbers. The operators (j)f and Ux have the 
action 



4)f\fi) = Uxlfi) = + A); 



(51) 



i.e., l/i) is an eigenstate of the smeared field operator 
and Ux is the generator of field translation. 

With this realization it is evident that configuration 
eigenstates are normalizable. This is one of the main 
difference between the polymer and Schrodinger quanti- 
zation schemes. It is because of this that the momentum 
operator does not exist in this quantization, but must be 
defined indirectly using the translation generators by the 
relation 



Quantum energy density 



To compute the expectation value in the semiclassi- 



cal constraint (45) we use the Gaussian coherent state 
peaked at the phase space values ((/)o,P0): 



IV') 



1 °° 



Af 



= ^-i<Pk-<Por/2cr^^-iP<p<PkVo^ 



(53) 



where = Afc/Fo^^^ is an eigenvalue of the scalar field 
operator derived from ^/ in Eq. (51). This state is cho- 
sen because it ensures that a large universe satisfies the 
classical Einstein equations, but introduces quantum cor- 
rections when the universe is sufficiently small [21]. 

Computing the expectation value of the operator p| in 
the above state we find 



{M\^) = 2A^ - exp(-eVs') cos(2e)] , (54) 



where 



e = XP^a- 



CTV0P4,, 



(55) 



are scale invariant variables; the first is just the expo- 
nent in the field translation operator, and the second is 
a measure of the with of the semi-classical state (since 
P(f) is a constant of also of the semiclassical equations of 
motion if the scalar field potential is zero). Hence the 
metric function on the boundary given by Eqn. (46) is 



m = 1 



^[l-exp(-eVS^)cos(2e)]„^^/^^. 

(56) 

It is straightforward to verify that in the limit A ^ 00 this 



(52) reduces to (36). The same is true for large R with fixed 



lambda. However for sufficiently small R the exponential 
factor goes to zero and f{R) is regular as ^ 0. This of 
course is just a consequence of the behaviour of the scale 
factor with the quantum matter source in the interior 
(which is communicated to the exterior via the junction 
conditions). It is this fact that significantly modifies the 
late stages of collapse in this model and gives singularity 
avoidance. 

A numerical integration of the star trajectory appears 
in Fig. 4 for the case k = 2. This is to be compared with 
Fig. 2 which is the same situation but with classical 
scalar field. Some salient features are that the junction 



8 



trajectory goes approaches R = zero exponentially indi- 
cating a rapid disappearance of mass rapid together with 
horizon formation and evaporation. Correspondingly, the 
metric function m{v) indicates a non-zero mass function 
that rises and falls to zero, and a variation of g{v) from 
negative to positive and eventually to zero. The weak 
energy condition is violated for the range of v values in 
which the apparent horizon shrinks 



trajectories have the same acceleration dpR/dv'^. This 
may provide an exactly solvable model where a star 
forms, becomes a black hole, and then evaporates. 

We thank Andreas Kreienbuehl and Sanjeev Seahra for 
discussions. This work supported by Natural Science and 
Engineering Research Council of Canada. 



CONCLUSIONS AND DISCUSSION 

We have constructed a new class of exactly solvable 
gravitational collapse models with two main results. The 
first is a classical generalization of the OS models to in- 
clude scalar field matter together with a clear interpre- 
tation of the exterior metric. A novel feature of this is 
that the exterior spacetime can be static or dynamic de- 
pending on the value of the equation of state parameter 
k in the stress-energy tensor. 

The second is a result in semiclassical gravity using 
the polymer quantization method, where we find that the 
behaviour of classical models is drastically modified due 
to matter quantization. Specifically for the k = 2 case 
we find that a singularity inside a horizon is replaced 
by formation and subsequent evaporation of an apparent 
horizon, with an exponentially disappearing mass. Al- 
though our model does not have all the features of a full 
field theoretical collapse model, such as the scalar field in 
spherical symmetry, it is interesting that it gives physi- 
cally desirable features of gravitational collapse that have 
been conjectured by various authors. 

The semiclassical model also provides a scenario to ex- 
plore other quantum states, values of /c, and the case 
of closed and hyperbolic interior metrics. An additional 
feature is that the junction equation (19) permits both 



expanding and contracting trajectories that permit a 
matching of these solutions. An example where this 
is possible arises in the closed universe case for which 
the parameter c in (19) must be less than one. This 



means that the trajectory reaches a turning point at a 
finite advanced time value (unlike the fiat case discussed 
here). This makes it possible to construct solutions that 
join collapsing and expanding branches at points where 
dR/dv = 0; at such points = (1,0,0,0) and the two 
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